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Abstract 

We confirm the Kosterlitz-Thouless scenario of the roughening transition 
for three different Solid- On- Solid models: the Discrete Gaussian model, the 
Absolute- Value-Solid-On-Solid model and the dual transform of the XY model 
with standard (cosine) action. The method is based on a matching of the 
*~ p 1 renormalization group flow of the candidate models with the flow of a bona fide 

KT model, the exactly solvable BCSOS model. The Monte Carlo simulations 
(— | are performed using efficient cluster algorithms. We obtain high precision 

estimates for the critical couplings and other non-universal quantities. For the 
XY model with cosine action our critical coupling estimate is /3^ Y = 1.1197(5). 
For the roughening coupling of the Discrete Gaussian and the Absolute-Value- 
Solid- On- Solid model we find K% G = 0.6645(6) and I(£ sos = 0.8061(3), 
respectively. 
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1 Introduction 

In 1951, Burton, Cabrera and Frank pointed out that a phase transition may occur in the 
equilibrium structure of crystal surfaces [1]. Such a phase transition from a smooth to a 
rough surface is called a roughening transition. 

The roughening transition was observed in experiments by Heyraud and Metois on 
metal crystallites with a diameter of only a few micrometres [2]. More recent results for 
the roughening transition of a (110) surface of Ag were reported by Robinson et al. [3]. 
They found a roughening temperature Tr = 790 ± 20 K. Several other groups have carried 
out beautiful experiments on roughening transitions for helium crystals in equilibrium 
with superfluid helium [4]. In these systems roughening temperatures of 0.35 K, 0.9 K 
and 1.28 K can be observed for the various inequivalent types of facets. Burton, Cabrera 
and Frank suggested that a growing layer of the crystal could approximately be described 
by the exactly soluble 2-dimensional Ising model. The idea was to represent the completed 
part of the layer by spin +1, while the vacancies are represented by spin —1. However, 
the model is quite unrealistic because it forbids the formation of a new layer before the 
previous one is completed. 

A more appropriate description of a crystal in equilibrium with its vapor is given by 
the 3-dimensional Ising model, where a spin +1 represents a site occupied by an atom, 
while spin —1 represents a vacancy. The boundary conditions are chosen such that an 
interface forms between a region where a majority of spins are +1 and a region with most 
spins equal to —1. This interface is the model of the crystal surface. In 1973, Weeks 
et al. performed a low-temperature expansion for the width of an (100) interface in a 
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3- dimensional Ising model on a simple cubic lattice with isotropic couplings [5]. They 
found a roughening temperature Tr = 0.57 T c . Here, T c denotes the critical temperature 
of the bulk phase transition. The approximation of the interface by the 2-dimensional 
Ising model yields Tr = 0.503 T c . 

A duality transformation exactly relates the 3-dimensional Ising model with the 3- 
dimensional Zi gauge theory. The roughening transition therefore also occurs in this gauge 
model [6, 7]. It was shown that also 4-dimensional lattice gauge models with continuous 
gauge group undergo a roughening transition [8-13]. The roughening transition is related 
to a restoration of symmetries broken by the introduction of the lattice. 

A fairly good approximation of the Ising interface is given by Solid- On- Solid (SOS) 
models. The idea is to neglect overhangs of the interface and bubbles in the bulk. The 
variables of the 2-dimensional SOS models represent the height of the interface with respect 
to some reference plane. 

SOS models are related with 2-dimensional XY models by a duality transformation [14]. 
The roughening transition of an SOS model corresponds to a transition of an XY model 
from a massive to a massless (spin wave) phase. In both representations, the transition 
is believed to be of the Kosterlitz-Thouless (KT) type [15]. The KT transition is a phase 
transition of infinite order, characterized by a very weak singularity in the free energy 
and an exponential singularity of the correlation length at the critical temperature. For 
reviews on the roughening transition see, e.g., [16]. 

Unfortunately for nearly all of the SOS models there is no rigorous proof that their 
phase transition is really of the KT type. For rigorous work on the existence of a phase 
transition from a massive to a massless phase, see [17]. See also [18], where the KT nature 
of the transition is put into question. 

In order to confirm or reject the KT nature of the roughening transition, many Monte 
Carlo simulations [19] were performed of various SOS models and of the interface in the 
3-dimensional Ising model [20-31]. 

Many Monte Carlo simulations were done to investigate the phase transition for the 
XY model with cosine or Villain action [32-49]. 

Most of the Monte Carlo studies of KT candidate models are based on a direct com- 
putation of critical quantities such as the correlation length or the susceptibility. 

By fitting the data with different ansatze one tries to rule out the power law singula- 
rities of conventional phase transitions. However, it turns out to be very demanding to 
get data for sufficiently large correlation lengths with good statistics. Only the results 
of the latest simulations using cluster algorithms really favour a KT transition against a 
second-order phase transition, while the estimate for the transition temperature still has 
a relative error of order 1% [50, 51]. 

There have also been attempts to study the KT scenario with the help of the Monte 
Carlo Renormalization Group (MCRG). See, e.g., [52]. 

In this paper, we attack the problem with a new method. 1 This is based on the fact 



J Part of this work is contained in [53]. A short report of our method and results is published in [54]. 
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that one of the SOS models, the BCSOS model, can be solved exactly [55-57]. The BCSOS 
model has been proved to exhibit a KT transition. The critical coupling is exactly known. 
In addition, the correlation length and other quantities can be computed exactly [55]. 

It was proposed long ago to improve the numerical study of SOS models by a comparison 
with BCSOS results [22]. In this report we give this comparison a precise meaning in the 
framework of the renormalization group (RG). We verify the KT scenario for several 
models - the ASOS model, the Discrete Gaussian (DG) model and the dual transform of 
the XY model with cosine action - by demonstrating that their long-distance RG flow at 
the critical point precisely matches with the flow of the critical BCSOS model. Stated 
differently, we demonstrate that the candidate models are in the same universality class 
(in the sense of Wilson's RG) as the BCSOS model. 

The matching is demonstrated by comparing Monte Carlo data for expectation values 
of "blocked correlation functions". All data are generated using cluster Monte Carlo 
algorithms that do not suffer from CSD or have strongly reduced CSD. References for 
cluster algorithms for spin models are [58] - [61]. [58-61]. Cluster algorithms for SOS 
models were introduced and studied in [62-64]. Our RG comparison is designed in such 
a way that finite-size effects are exactly cancelled. By simulations on reasonably small 
lattices we obtain results for the critical couplings, which are competitive in precision 
with estimates from much more expensive Monte Carlo studies. We also get estimates 
for the non-universal constants determining the asymptotic behaviour of the correlation 
length. 

The perhaps most important result of our study is, however, the demonstration that 
the models (with a high level of confidence) are in the same universality class as the BCSOS 
model. We consider this as an unambiguous confirmation that their phase transition is of 
the KT type. 

This article is organized as follows: In section 2 we introduce the SOS models and their 
most important properties. In section 3 we summarize the heuristic picture of the KT 
phase transition using the flow diagram of Kosterlitz and Thouless. In section 4 we define 
an RG transformation for theories on finite lattices. The RG flow defined in section 4 lies 
at the heart of the matching analysis presented in section 5. There, we demonstrate that 
all the considered SOS models share the same universality class with the BCSOS model. 
Appendix 1 is devoted to the solution of the KT equations. In Appendix 2 we present the 
exact computation of the finite-size RG for the Gaussian model. 

2 Solid-on-solid models 

Solid-on-solid (SOS) models are of interest both theoretically and in the study of crystal 
interfaces. In this section we define some of these models and give an account of their 
most important properties. For reviews on SOS type of models we refer to [16]. 

All SOS models have in common that they are 2-dimensional lattice spin models. The 
spins h x take values in an unbounded discrete set S (isomorphic to the integer numbers). 
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The interaction energy (Hamiltonian) is invariant under global shifts h x — >■ h x + M for 
M G S. A typical partition function for such a model looks like 2 

Z = ]Texp I- £ V(h x -h y )\ . (1) 

{/i} \ <x,y> j 

In this example, the Hamiltonian is a sum of contributions depending on pairs of nearest- 
neighbour spins only. The summation in eq. (1) is over equivalence classes of spin confi- 
gurations {h}. The classes are defined by identifying two configurations that differ only 
by a global shift M G S. 

Our first example of an SOS model is the DGmodel: It is of the type defined in eq. (1) 
with 

V DG = K DG (h x - hy f . (2) 

The spin variables h x take integer values. Note that the Hamiltonian looks exactly like 
that of a continuous Gaussian model. However, the restriction of the h x to integer values 
introduces a non-trivial interaction. Let us interpret the h x as heights with respect to a 
certain base. For finite positive K DG the Hamiltonian will favour that neighbouring spins 
take similar values. When K DG is large enough, the surface will not fluctuate too wildly. 
As a consequence one expects that the surface thickness squared, 

a 2 = lim ((h x -h y ) 2 ), (3) 

\x — y \— >oo 

is finite: the system is in the "smooth" phase. On the other hand, if K DG is below a 
certain critical value, the surface becomes "rough", and the surface thickness diverges. 

In [17] it was proved for a class of SOS models that the 2-point correlation function 
((ho — h x ) 2 ), for sufficiently small coupling K , goes like ln(|a;|) at large distance x. 

Furthermore, it follows from convergent cluster expansions that ((ho — h x ) 2 ) stays 
bounded for all x if K is sufficiently large. 

On finite lattices with N = L X L sites we define 

ct2 = ^2E<(^-M 2 >- (4) 

The surface thickness defined through this equation behaves exactly like the surface thick- 
ness defined via eq. (3): It stays finite for L — > oo if K DG is above a critical value, and it 
diverges with increasing L if K DG is below this critical value. In this work, we shall refer 
to definition (4) throughout. The transition between the two phases is called roughening 
transition. The theory of Kosterlitz and Thouless makes detailed predictions about the 
nature of this transition, see section 3. 



2 A factor 1/ksT, where ks denotes Boltzmann's constant and T the temperature, is absorbed in the 
definition of the interaction. 
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The Discrete Gaussian model is dual to the XY model with Villain action [14]. This 
model is defined by the partition function 

Z v = r ]JdQ x n B(Q X -Q y ), (5) 

" v x <x,y> 

with 

oo 

B(&) = J2 exp(-i/3 v (0-27rp) 2 ) (6) 



p= — CO 



and 



1 K DG . (7) 



2(i v 

The index "V" here refers to "Villain". 

The XY model with "standard (cosine) action" has the partition function 



zxy = r n rf0 - ex p [p xy e 

n x y <x,y> 



cos(Q x - Q y ) 



The standard action is the action discussed mostly for an XY model. The dual of this 
model is given by the partition function 

{h} <x,y> 

where the I n are modified Bessel functions. Again h x is integer. 

The Absolute- Value-Solid-On-Solid (ASOS) model is the SOS approximation of an in- 
terface in an Ising model on a simple cubic lattice on a (OOl)-lattice plane. It is defined 

by 

Vasos = K ASOS \h x - h y \ . (10) 

We shall finally introduce the BCSOS (Body Centered Solid-On-Solid) model. It will 
play a prominent role in this paper. The BCSOS model was introduced by van Beijeren 
[65] as an SOS approximation of an interface in an Ising model on a body-centered cubic 
lattice on a (OOl)-lattice plane. The effective 2-dimensional lattice splits in two sublattices 
like a checker board. In the original formulation, on one of the sublattices the spins take 
integer values, whereas the spins on the other sublattice take half-integer values. We adopt 
a different convention: spins on "odd" lattice sites take values of the form 2n + |, and 
spins on "even" sites are of the form 2n — |, n integer. As a consequence, the effective 
distribution for block spins (defined as averages over square blocks of size B X B, where 
B is an integer) will be centered around integer values instead of half-integer values, and 
the configurations with minimal energy will have integer average values, as is the case for 



SOS models 



7 



K 


^asympt 


€exact 


K 


^asympt 


€exact 


1.00 


1.577 


2.033 


0.20 
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0.3080- 10 9 


0.3156- 10 9 



Table 1: Comparison of ^ BCSOS from the asymptotic and the exact formula 



the other SOS models defined in this section. The partition function of the BCSOS model 
can be expressed as 

Zbcsos = £exp ( -K BCSOS J2 \K ~ K\ ) . (11) 

{h\ \ [x,y] / 

The sum is over next-to-nearest-neighbour pairs [x, y], and nearest-neighbour spins h x and 
h y obey the constraint \h x — h y \ = 1. In [65] Van Beijeren showed that the BCSOS model is 
isomorphic to the F-model, which is a special six- vertex model. The configurations of the 
BCSOS model are in one-to-one correspondence to the configurations of the F-model. The 
F-model can be solved exactly with transfer matrix methods [56, 57, 55]. The roughening 
transition occurs at 

k bcsos = l m2 . (12) 
The exact formula for the correlation length is [55]: 

1 J 1/2 ~ / l + x ^ \ 2 \ 

ZBCSOS - " ln j 2 * ^1 ^ + ^-2 J j . 

x = exp (— arcosh (I exp(41i') — l)) . (13) 

For K \ Kr, the correlation length behaves like 

t BCSOS ^\e*J-^ K -h) , (14) 

The essential singularity of the correlation length at the critical coupling Kr is typical 
for a KT model. It is instructive to compare the results of the exact and the asymptotic 
expression in the neighbourhood of the critical point. Table 1 shows that at correlation 
length 189 there is still a relative deviation of 7 % between the exact result and the result 
from the asymptotic formula. This nicely explains why the determination of amplitudes 
(occurring as parameters in the asymptotic correlation length formula) from correlation 
length measurements in KT models is so difficult. 

The free energy per volume is also exactly known. For K > Kr, 



f bcsos = -l n Z BCSOS /volume 
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2K [i A , exp(-mA)sinh(mA) | 
\ 2 ,~ L m cosh(mA) J ' 
A = arcosh (± exp(4ii.') - l) . (15) 

For K < JtfR, the free energy f BCSOS has the following integral representation: 

rBCSOS r> ts 1 f°° dx n ^cosh a; — cos 2/i 



Jo cosli[irx /2/j,) \ cosh a; — 1 

(16) 

fj, = arccos (| exp(41i') — l) . (17) 

Figure 1 shows f BCSOS ' m the neighbourhood of the roughening coupling Kr = |ln2. 
The free energy and all its derivatives stay finite at the transition point (when suitably 
defined as limits from the right and from the left, respectively). However, there is an 
essential singularity: This can already be seen if one directly inserts K = Kr in eqs. (15) 
and (16). In both cases one obtains ln2, to be compared with the right-left limit in fig. 
1. It can be shown [55] that the singular part of the free energy in the vicinity of the 
roughening point is of the form 

f s lf OS - exp ' « = ^ ■ (18) 



3 Kosterlitz-Thouless theory 

The Kosterlitz-Thouless theory is conveniently discussed in the framework of the Sine 
Gordon model, a system with continuous spins and lattice Hamiltonian 

n ^ = Yd ^ (^-^) 2 -^E cos ( 27r ^)- (19) 

' <x,y> x 

The coupling constant z is called "fugacity", because the model can be exactly transformed 
into a 2-dimensional Coulomb gas, where z plays the role of the fugacity for charged 
particles. In the limit of vanishing fugacity the roughening coupling of the Sine Gordon 
model is 

lim/3f^) = -. (20) 

z— >0 7T 

The critical behaviour of the model is determined by the long-distance properties of 
the correlation functions. The most systematic way to analyze long-distance properties 
of a statistical model is the renormalization group as pioneered by Kadanoff et al. [66] 
and worked out systematically by Wilson [67]. At the heart of the renormalization group 
approach to critical phenomena lies the notion of the effective Hamiltonian, which encodes 
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Figure 1: The free energy per volume of the BCSOS model as a function of 



the properties of the model at length scales larger than the (microscopic) scale of the 
fundamental Hamiltonian. For a lattice model the effective Hamiltonian 7i e ff is defined 
through 3 



The effective theory lives on a lattice A', which has lattice spacing B (in units of the 
original lattice spacing). The block spins (f> x i are defined as averages of the original spins 
(p x of square blocks x' of size B X B. 

In general, one needs a non-trivial rescaling of the block spin field (in accordance 
with a non- vanishing critical exponent r]). However, for the Z-symmetric models studied 
in this work, there is no anomalous dimension: Note that for SOS models the effective 
Hamiltonians defined through eq. (21) also have a global Z-symmetry. This property 
would be destroyed by introducing additional rescaling in the definition of block spin 
variables 4 . 

A renormalization group flow is defined by iterating the renormalization group transfor- 
mation. In the case of a linear block spin definition the iteration of a single step is exactly 

3 There exist, of course, many other possible definitions of an effective Hamiltonian. 

4 The critical fixed point, which is a continuous Gaussian model, has even larger symmetry, namely the 
real numbers R. This phenomenon is called symmetry enhancement. 




(21) 
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equivalent to going to larger and larger block size. The property R{B\) i?(i?2) = R{B\ B2) 
is called half group property. 

One of the practical difficulties one has to deal with is the proliferation of coupling 
constants: the exact effective Hamiltonian can no longer be exactly parametrized by a 
finite number of coupling constants. There is an infinite number of effective coupling 
constants. However, one assumes that only a few couplings are important for the critical 
behaviour. This has of course to be demonstrated. The study of the effects from truncating 
the effective Hamiltonian is a difficult problem. 

KT theory is based on the assumption that the RG flow of the critical or nearly critical 
Sine Gordon model (and also of other KT models) at long distance is well described by 
two parameters (3 and z [15]. 

The study of the flow of (3 and z with increasing length scale is most convenient in a 
continuum formulation, where the lattice cutoff is replaced by a momentum space cutoff. 
A length scale is then identified with the inverse momentum cutoff. The cutoff can be 
changed continuously. One can derive differential equations for the flow of the effective 
Hamiltonian [15]. The KT flow equations are 

y = -xy , 

x = -y 2 , (22) 

where the dot denotes the derivative with respect to the scale parameter t. The corres- 
ponding length scale is B = exp(i), and x and y are related to (3 and z by 

x = ir(3 — 2 , 

y = cz, (23) 

where the constant c is non-universal, i.e. depends on the choice of the regularization 
procedure. The equations can be exactly solved, see Appendix 1. The solutions satisfy 

the equation y 2 = E + x 2 , or 

y = \/E + X 2 . (24) 

Typical renormalization group trajectories are shown in fig. 2. The flow diagram encodes 
the full information about the critical behaviour. Region I (E < 0) corresponds to the 
massless phase. Here, the trajectories run exponentially fast (in the variable / = In B) to 
y = 0: the long-distance behaviour is that of a massless Gaussian model. The point where 
a trajectory hits the s-axis corresponds to a /3-value that is called /5 e ff . It is the /3-coupling 
of the free field theory describing the long distance behaviour of all models that "come 
down" the trajectory. Note that this is a whole class of models, which differ only in the 
"time" they need to reach the s-axis. 

Regions II (E > 0) and III (E < 0) correspond to the massive phase. Here, with 
increasing length scale, the fugacity increases, and the global symmetry under shifts ip x — f 
(f x + integer is spontaneously broken. The separatrix between regions I and II is the critical 
line. Models on this line are driven into the fixed point (x,y) = (0,0), which corresponds 
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KT renormalization group trajectories 
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Figure 2: Renormalization group trajectories in a Kosterlitz-Thouless model. Region I 
(E < 0) corresponds to the massless phase. Here, the trajectories run exponentially fast (in 
the variable t) to y = 0: the long-distance behaviour is that of a Gaussian model. Regions 
II (E > 0) and III (E < 0) correspond to the massive phase. Here, with increasing length 
scale, the fugacity increases, and the global symmetry under shifts ip x — f ip x + integer is 
spontaneously broken. The separatrix between regions I and II (E = 0) is the critical line. 
Models on this line are driven into the fixed point (x,y) = (0,0) which corresponds to 
fi = 2/-7T and z = 0. The meaning of the dotted line intersecting the trajectories in regions 
I and II is explained in the text. 
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to (3 = 2/-7T and z = 0. Note that the approach to the fixed point goes like ~ t~ x on the 
critical trajectory. 

From the flow diagram one can derive (see, e.g., [68]) that the correlation length of a 
KT model diverges like 



Here, A, C, and (3 C are non-universal constants, i.e. they have different values for different 
KT models. Compare with eq. (14). 

The KT equations for the Sine Gordon model are derived using several approximations; 
the most important of these is that the fugacity is small. However, if one assumes that 
(x, y) = (0, 0) is the relevant fixed point for the critical behaviour, then the trajectories of 
a critical or nearly critical KT model will eventually enter the region where the fugacity 
is small and where the KT equations become valid. Depending on the details of the 
microscopic Hamiltonian, the "time" before the effective Hamiltonian comes close to the 
fixed point can be very long, i.e. the flow close to the critical line is very slow. This is one 
of the reasons why the unambiguous confirmation of the KT scenario for realistic models 
is so difficult. 

Idealizing things very much, one of the SOS models introduced in section 2 can be 
represented by a line like the dotted line in fig. 2 intersecting the trajectories in regions I 
and II. The line should be carefully interpreted as follows: to each value of the SOS coupling 
K , there corresponds one of the trajectories in the flow diagram the model eventually will 
run on after several RG steps. 

4 Finite lattice renormalization group 

Universal properties of a statistical system do not depend on short distance details, but 
only on the nature of long wavelength fluctuations. This suggests to remove the irrelevant 
high frequency degrees of freedom by applying a coarse graining (block spin) procedure as 
introduced in section 3. 

Universality, which was first introduced as the coincidence of the critical indices of 
various models, can be expressed as a convergence of the renormalization group flow to a 
universal flow as T — > Tr and the number of block spin transformations goes to infinity. 
Here we use T as a representative for any coupling that is driven to a critical value in 
order to make the correlation length diverge and the system become critical. 

Consider two models with different microscopic Hamiltonians that belong to the same 
universality class, i.e. that have the same critical indices. The above statement says that 
if the two models are both at criticality then their effective Hamiltonians will converge 
towards the same fixed point Hamiltonian. It might even happen that the two flows of 
Hamiltonians will come close to each other already a long time before they are really close 
to the fixed point. (We shall actually observe this phenomenon for the models studied in 




(25) 
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this paper.) Note, however, that the two systems might need a different number of RG 
steps to reach a certain point on the universal trajectory. 

It is suggestive to compare renormalization group flows in order to test for universality 
properties. To directly compare effective Hamiltonians, one would have to parametrize 
them in terms of coupling constants. However, there are infinitely many couplings already 
after a single renormalization group step. It turns out to be very difficult to determine the 
coupling constants of the blocked system by analytical calculations [69] or Monte Carlo 
simulations (MCRG) [70]. 

We shall deal with the problem in a similar fashion as Shenker and Tobochnik [71] 

did for the 2-dimensional 0(3) model and Wilson [72] for the 4-dimensional SU(2) gauge 

model. We define effective Hamiltonians for finite lattices A that consist of / X / blocks, 

(l B) 

where each of the blocks contains B X B sites. The effective Hamiltonian % ft defined 
through 

e X p(-^ ) (^>)) = / n d<p x ^ V {-H{<p)) I] 5 \<t>x- - B- 2 <Px\ ( 26 ) 
xEA x'£A' \ xEx' j 



is a function of I 2 block spin variables <f> x i. A renormalization group flow is now defined 

(l B) 

as the sequence of effective Hamiltonians H ft , for fixed / and increasing B. Note that 
different Z's lead to different flows. 

(/ B) 

In order to monitor the flow of the % ft , we do not necessarily have to compute the 

(l B) 

still infinitely many couplings in % ft ■ We can instead consider a set of suitably chosen 
observables, 



E {i,B) 



(M<t>))l,B = Z- 1 j D4> exp(-?4ff B) ) M<f>) . (27) 

(/ B) 

A convergence of the flow of effective Hamiltonians % ft towards a fixed point will imply 

(/ B) (I B) 

also the convergence of the flows of the E\' . The crucial point is that the E\ ' can 

be expressed as expectation values in the original system with Hamiltonian H on a lattice 

with size L = IB. One just has to measure correlation functions of block averages of the 

original system, that can be simulated, e.g. with efficient Monte Carlo algorithms. 

(I B) 

It is worth noting that the couplings E\ ' introduced above can be interpreted as 

phenomenological couplings as introduced by Nightingale and by Binder, see [73]. 

(I B) 

We conclude this section by presenting results for the flow of the % ft for the free 
massless field theory (Gaussian model) in two dimensions. This may serve as an illustration 
for the convergence of the flow towards a fixed point. The Hamiltonian is 

H(<P) = \{<P, "A?) =\ E (<P*~ Vy? ■ (28) 

<x,y> 

For this theory the effective Hamiltonian can be computed exactly (see Appendix 2), 

*2r fl) (*) = \ (*> " V) = \ E [" A2r fl) ] x , y , <t>y • (29) 

x' ,y' 
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-0.27744 


-1.49801 


-0.03343 


0.06616 


128 


-10.1783 


3.58908 


-0.27815 


-1.49950 


-0.03315 


0.06613 


256 


-10.1793 


3.58964 


-0.27833 


-1.49988 


-0.03308 


0.06612 


512 


-10.1796 


3.58978 


-0.27837 


-1.49997 


-0.03307 


0.06612 


1024 


-10.1796 


3.58982 


-0.27838 


-1.49999 


-0.03306 


0.06612 


2048 


-10.1797 


3.58983 


-0.27838 


-1.50000 


-0.03306 


0.06612 



Table 2: Flow of the effective Laplacian for / = 4. The components are 
arranged as indicated in eq. (30) 

(l B) 

Table 2 shows the components of 1 on a 4 X 4 lattice, which are arranged according 
to the following scheme: 

a 
b 
d 
b 

5 RG matching with the BCSOS model 

Using the cluster algorithm described in [63], we simulated the BCSOS model at the rou- 
ghening coupling Kft CSOS = | In 2 on square lattices of size L X L with periodic boundary 
conditions. For a list of the used i's and the statistics, see table 3. There we give as an 
example the results for the squared interface width a 2 defined in eq. (4). 



L 


a 2 


stat 


L 


a 2 


stat 


12 


1.02860(32) 


2.8 


48 


1.34674(34) 


2.2 


16 


1.09624(42) 


1.4 


64 


1.41080(26) 


4.4 


24 


1.18996(40) 


1.6 


96 


1.50038(30) 


5.3 


32 


1.25464(38) 


1.5 


128 


1.56362(32) 


4.9 



Table 3: Results for the squared interface width a 2 as function of the lattice 
size L for the critical BCSOS model. The statistics stat is given in units of 
10 6 single cluster updates 

Figure 3 shows a 2 — (2/7r 2 )ln L as a function of InL. KT theory predicts that this 
quantity should converge towards a constant for large L, see [15]. The figure shows a 



b 


d 


b 


c 


e 


c 


e 


f 


e 


c 


e 


c 



(30) 
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0.52 



2.0 2.5 3.0 3.5 4.0 4.5 5.0 

InL 



Figure 3: The squared surface width a 2 of the BCSOS model minus its large 
L behaviour as anticipated from KT theory 



significant deviation from this behaviour. This means that for the lengths fitting on 
lattices with L < 128 the effective fugacity is not small. This is a consequence of the fact 
that along the critical line the flow towards the Gaussian fixed point is very slow (like 
~ 1/lnL, see section 3). 

Our method to monitor the RG flow is to compute the flow of blocked observables (see 
section 4). The lattice is divided into / X / blocks of size B X B, with / = 1,2, 4. 5 Linear 
block spins (f> x i are defined according to eq. (21). 

Motivated by KT theory we measured two types of block observables: those that are 
sensitive to the flow of the kinetic term (flow of K), and those that are sensitive to the 
fugacity. For the first type of observables we chose 

A ^ = (^p £ <t>v>) 2 ) , (31) 

\ <x',y'> i 
where < x', y' > are nearest-neighbour pairs on the block lattice, and 




5 A posteriori we found that measuring also I > 4 data would have been useful. However, the small 
l-values enabled us to save all block spin configurations on disk. This gave us flexibility in the data analysis. 
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L 


,(0) 
^1,2 


,(0) 
^2,2 






8 


n 

w 


1 3671 Q 


n 

w 




0.293527 


0.380581 


12 


n 

w 


1 9R791 


n 

w 




0.257225 


0.340481 


16 


n 

w 


1 931 47 


n 

w 


1 71 87^ 

1 1 lulu 


0.243918 


0.326090 


24 


n 

U 


1 90565 


n 

U 


1 68981 


0.234147 


0.315663 


32 


n 

U 


1 1 9655 


n 

U 


1 67969 


0.230662 


0.311978 


48 


n 


1 1 9009 


n 


1 67945 


0.228148 


0.309333 


64 


o 


118773 


o 


166992 


0.227263 


0.308404 


96 





118609 





166811 


0.226628 


0.307739 


128 





118551 





166748 


0.226406 


0.307507 


256 





118496 





166687 


0.226191 


0.307282 


512 





118482 





166672 


0.226137 


0.307226 


1024 





118479 





166668 


0.226124 


0.307212 


2048 





118478 





166667 


0.226120 


0.307208 


4096 





118478 





166667 


0.226119 


0.307207 


8192 





118478 





166667 


0.226119 


0.307207 



Table 4: Exact results for A\ I and A 2 



where [x', y'] are next-to-nearest-neighbour pairs. These quantities are defined for / > 1. 
For the actual matching procedure to be described below, we also employed the quantities 

A {0) \ B = 

Ay = (0) °° Ai t , for / = 1,2. (33) 

A, i 

Here, the A^ denote the same quantities as defined in eqs. (31) and (32), taken, however, 

in the continuous Gaussian model with the Hamiltonian defined in eq. (28). The A^ can 

be computed exactly, see Appendix 2. In table 4 we give the values of A^j and A^j 
for L < 8192. Within the accuracy obtained (6 digits), the infinite B limit is reached for 
L = 4096. 

As a monitor for the fugacity we chose the following set of quantities (defined for 
/= 1,2,4): 

-43,/ = ^^cos(l -271-^/)^ , 

Aa,1 = (j^ ^COs(2 • 27T^/)^ , 

As,, = ^X>s(3-27r<^ • (34) 

We believe that all important information about the large- distance RG flow of an SOS 
model close to or at criticality can be monitored by these observables or by a subset of 
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L 


/ 


5 


A i,i 


A 2 ,i 


A 3J 


A 4 ,i 


A 5 ,i 


12 


1 


12 






0.2842(17) 


0.0777(14) 


0.0229(13) 


16 


1 


16 






0.2655(27) 


0.0669(20) 


0.0212(18) 


24 


1 


24 






0.2383(31) 


0.0551(22) 


0.0149(19) 


32 


1 


32 






0.2261(36) 


0.0511(24) 


0.0105(20) 


48 


1 


48 






0.2056(36) 


0.0394(23) 


0.0077(18) 


64 


1 


64 






0.1969(30) 


0.0372(18) 


0.0051(13) 


96 


1 


96 






0.1821(37) 


0.0296(20) 


0.0052(14) 


128 


1 


128 






0.1719(43) 


0.0283(21) 


0.0058(15) 


oo 


1 


oo 






0.0 


0.0 


0.0 


12 


2 


6 


0.08891(16) 


0.11983(26) 


0.2403(10) 


0.0668(7) 


0.0227(6) 


16 


2 


8 


0.08569(22) 


0.11745(36) 


0.2258(16) 


0.0582(10) 


0.0157(9) 


24 


2 


12 


0.08319(22) 


0.11505(37) 


0.2009(19) 


0.0445(11) 


0.0101(9) 


32 


2 


16 


0.08187(23) 


0.11395(39) 


0.1892(22) 


0.0390(11) 


0.0084(10) 


48 


2 


24 


0.08114(22) 


0.11343(37) 


0.1746(22) 


0.0328(11) 


0.0058(9) 


64 


2 


64 


0.08097(17) 


0.11339(29) 


0.1655(18) 


0.0263(8) 


0.0042(7) 


96 


2 


48 


0.08076(20) 


0.11326(33) 


0.1528(22) 


0.0239(9) 


0.0046(7) 


128 


2 


64 


0.08038(21) 


0.11268(35) 


0.1435(26) 


0.0208(10) 


0.0033(7) 


oo 


2 


oo 


0.075425 


0.106104 


0.0 


0.0 


0.0 


12 


4 


3 


0.19524(17) 


0.24472(24) 


0.2438(5) 


0.1121(4) 


0.0586(4) 


16 


4 


4 


0.17686(21) 


0.23046(31) 


0.2226(7) 


0.0781(5) 


0.0407(5) 


24 


4 


6 


0.16521(20) 


0.21978(30) 


0.1933(8) 


0.0535(6) 


0.0177(5) 


32 


4 


8 


0.16087(21) 


0.21539(31) 


0.1793(9) 


0.0438(6) 


0.0122(5) 


48 


4 


12 


0.15798(18) 


0.21278(28) 


0.1632(9) 


0.0349(5) 


0.0087(4) 


64 


4 


16 


0.15646(14) 


0.21136(21) 


0.1522(8) 


0.0296(4) 


0.0064(3) 


96 


4 


24 


0.15522(16) 


0.20995(25) 


0.1395(9) 


0.0248(5) 


0.0049(4) 


128 


4 


32 


0.15471(17) 


0.20928(27) 


0.1322(11) 


0.0219(5) 


0.0045(4) 


oo 


4 


oo 


0.143952 


0.195574 


0.0 


0.0 


0.0 



Table 5: Finite lattice renormalization group flow of the Aij for the critical 
BCSOS model 



them. The result for the flow of the A's for the critical BCSOS model is summarized in 
table 5. We also give the exact limits that these quantities should approach when the 
block size B is scaled to oo. KT theory predicts that A^^, A^\ and A$ t i have to converge 
to zero. The quantities A\ t i and A^^i are predicted to converge to 2/ it times the B — > oo 
limit of the same observables in the free field theory. These limits can be read off from 
table 4. A close look at the data reveals that even for large B the A's are still off their 
fixed points values. However, we shall see in the following that it does not matter that in 
the flow of the BCSOS data the fixed point is still somewhat away: The RG matching will 
take place a long time before the fixed point is close. Irrelevant couplings die out with a 
power of the length scale. The flow of the couplings is rapidly reduced to a 1-dimensional 
manifold. Along this remaining line, the fugacity dies out logarithmically with the length 
scale. Eventually the Gaussian fixed point with zero fugacity is reached. Since we know 
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the 1-dimensional manifold from the BCSOS model, we just have to recover it in the other 
models. This can be done even far away from the fixed point. 

The simulations of the DG model, the ASOS model and the dual of the XY model were 
also performed on quadratic lattices with periodic boundary conditions. We used the very 
efficient VMR cluster algorithms [62]. The same blocking prescription as for the BCSOS 
model was employed. The lattice sizes and couplings involved will be specified below. 

5.1 Determination of the roughening couplings 

There are two parameters that have to be tuned in order to match the RG flow of one of 
the SOS models with that of the critical BCSOS model. One can vary the coupling K sos 
of the SOS model. This allows one to walk on the approximate starting line in fig. 2. The 
flow of the SOS model can only match that of the critical BCSOS model if K — KpP^ . 
On the other hand, one can vary the ratio of the lattice sizes of the SOS model and the 
BCSOS model and, as a consequence, the ratio of the block sizes = B sos / B BGSOS . 

A ratio b^ s ^ 1 turns out to be necessary to compensate for the different positions of 
the approximate starting lines in fig. 4. 

Before we turn to the question of how to determine Kf> OS and b^ s in practice, let us 
write down the general condition for one of the SOS models to be in the same universality 
class as the critical BCSOS model. 

Matching condition: Universality holds if there exists a b m and a Kf> OS such that for all 
Af? s (b s m os B BCSOS , K s ° s ) = A BCSOS (b bcsos , K bcsos ) (35) 



in the limit of large B BGSOS , anc [ the corrections are of order (B BCSOS )~ UJ , with u > 
the leading correction to scaling exponent. As we shall see below, moderate B BGSOS are 
sufficient in practice. 

In order to determine the SOS roughening coupling Kff S from the RG flow data we 
proceeded as follows: For fixed values of L sos and for fixed /, we considered the two 
equations: 

j\SOS ^b sos K sos ^j = a bgsos ^b bcsos K BGSOS 
Al° s (B sos Jii os ) = A B f sos (B BCSOS ,K BCSOS ) . (36) 

We chose A\ and A3 because we consider these to be the most important observables 
for the monitoring of the RG flow. For each of the available values of B BGSOS listed in 
table 5 we solved these two equations for the couplings Kf os . To be able to do this we 
needed the A's for a range of couplings. We simulated the SOS models at the (to that 
time) best known estimate for their roughening coupling. The expectation values in a 
neighbourhood of the simulation point could then be obtained by extrapolating using a 
reweighting method [74]. 
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Figure 4: Matching of the finite size RG flow of two different SOS models in 
the KT flow diagram. Two successive points along the discrete trajectories 
are separated by a fixed scale factor. (In reality the models do not 'start' 
in the KT diagram but rather in a higher dimensional space of coupling 
constants.) Note that one of the models is one step 'ahead' of the other one. 
This explains that one needs the offset factor b m in order to match the flows. 
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Figure 5: Determination of f3 x Y for L XY = 32, L BCSOS = 32 and / = 2. The 
solid curve gives A XY as a function of /3 XY . The dashed curves indicate the 
statistical error. The solid straight horizontal line gives A BGSOS at K BCSOS . 
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For the determination of , we did not use A\ t i directly, but the 'improved' quan- 

tity Di t f. matching of the A^-flows of two models happens if and only if also the Dij- 
flows match. This is so because for B — > oo the factors A^(B) converge to fixed points 

A\ 0) (B = oo). This proves that we are allowed to use the Di i instead of the An without 
losing any control on the RG flow. Furthermore, the Dij flows converge more rapidly than 

those of the A^;. This will be demonstrated below. Roughly speaking, the A$ factors 
cancel irrelevant terms in the flow that anyway die out under successive RG steps. These 
terms are strong as long as B is small, and are partly due to discretization details, e.g. of 
the lattice Laplacian. We want to stress the point that the use of the _D's instead of the 
A's is by no means necessary: the changes on the larger lattices are negligible. However, 
using the D's allows one to observe a collapse to a universal trajectory on much smaller 
lattices. 

The solution of eq. (36) yields, for each (L,l) pair, two values: Kf os and K^ os . For 
an illustration of this first step, see fig. 5. Note that Kf os and K^ os will in general not 
be identical: one can expect matching only for a specific ratio 6^ os '. 

In a second step we plotted the values of Kf os and K^ os as a function of B BGSOS . 
The couplings were linearly interpolated in log B BC sos . The intersection of the two 
curves xf os (B BCSOS ) and K^ os (B BGSOS ) then uniquely determines an estimate for 
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Figure 6: Determination of (3p[ for L = 32 and / = 2. The curves 
give /3i Y and fif Y as functions of Bbcsos- The intersection of the curves 
uniquely determines f3^ Y and the matching Bbcsos- 

the roughening coupling Kf> OS of the SOS model, see fig. 6. 



matching in the sense described above. The results for the matching for the three models 
are summarized in table 6. 

For the three models, the results for the roughening coupling Kr obtained for the 
various lattice sizes L and sizes / of the blocked system are consistent with each other 
within statistical errors. Only the couplings for / = 4 on the smallest two lattice sizes and 
for / = 2 on the smallest lattice size deviate slightly from the rest. For the ratio of the 
matching block sizes b m = B sos / B BCSOS the observation is the same. This indicates an 
extremely fast convergence to a universal RG flow of the models, since even for such small 
block sizes as B BGSOS = 16 no deviation from the universal flow can be observed within 
our quite good statistics. To give estimates for the roughening coupling Kr for the three 
models we averaged the values obtained for the largest L with / = 2 and / = 4, and the 
second largest L only with / = 2. We arrive at the following results 





1.1197(5), 
0.6645(6), 
0.8061(3). 



K 



,-ASOS 



(37) 



The quoted errors are statistical, but according to the discussion above, the systematic 
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Dual of XY model 


L 


Pr, I = 2 


/?*,/ = 4 




b m , I = 4 


16 


1.1220(12) 


1.1257(8) 


0.84(5) 


0.75(1) 


24 


1.1214(13) 


1.1225(8) 


0.91(7) 


0.84(2) 


32 


1.1211(12) 


1.1214(8) 


0.93(10) 


0.85(3) 


48 


1.1199(11) 


1.1205(7) 


0.89(13) 


0.89(3) 


64 


1.1212(11) 


1.1201(7) 


0.89(5) 


0.82(12) 


96 


1.1189(11) 


1.1194(7) 


0.89(12) 


0.95(7) 


DG model 


L 


/?*,/ = 2 


/?*,/ = 4 


i ^ — 2 


b m , I = 4 


12 


0.6627(16) 


0.6607(13) 


0.31(4) 


0.40(2) 


16 


0.6650(13) 


0.6632(10) 


0.32(5) 


0.34(2) 


24 


0.6633(16) 


0.6645(8) 


0.30(6) 


0.34(2) 


32 


0.6650(16) 


0.6647(8) 


0.28(5) 


0.32(2) 


ASOS model 


L 


/?*,/ = 2 


/?*,/ = 4 


i ^ — 2 


b m , I = 4 


32 


0.8052(4) 




2.3(2) 




64 


0.8061(6) 


0.8058(3) 


2.8(3) 


2.4(1) 


128 


0.8061(6) 


0.8060(3) 


2.7(6) 


2.6(2) 


256 


0.8060(5) 


0.8062(3) 


2.8(6) 


2.9(4) 



Table 6: K R and b m = B sos / B BCSOS for the three SOS models as obtained 
from the matching of A\ t i and ^3^. 
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ones due to deviations from the universal parameter flow, should be much smaller. For 
the b m we find in a similar way 



b XY 


= 0.89(5), 




b DG 


= 0.31(2), 




t ASOS 

m 


= 2.8(3). 


(38 



5.2 Demonstration that the matching is unique for all Aij 

We want to demonstrate that all quantities Aij (Di,i) converge towards a universal flow 
with increasing B sos , provided that the couplings K sos are tuned to their critical values 
quoted in eq. (37), and that the block size ratios b^ s = B sos / B BCSOS are taken to be 
the matching values quoted in eq. (38). Recall that the critical couplings and 6 m 's were 
determined by imposing the matching condition for D\ t i and A^^ alone. 

The first task was to evaluate the observables A2J, A^\ and A$ t i at the critical cou- 
plings Kr determined above using D\ t i and only. The results are summarized in 
tables 7 to 9. For the sake of completeness we also give the values of the Dij, in table 
10. We furthermore demonstrate the universal matching by plotting all measured block 
observables at criticality as functions of the matching block size B BGSOS = B sos /b!^ s . 
The reader is invited to look carefully at figs. 7 to 23. To correctly interpret the plots, it 
is necessary to realize that the scale of the y-axis might differ from plot to plot. See, for 
example, figs. 7 and 8: The collapse of D\p onto a universal curve is much faster than for 
the corresponding quantity A\^- 
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L 


/ 




A 2 ,i 


A 3J 


A 4J 


A 5 ,i 


16 


1 






0.2601(19) 


0.0669(17) 


0.0175(13) 


24 


1 






0.2381(28) 


0.0497(20) 


0.0110(17) 


32 


1 






0.2260(22) 


0.0493(18) 


0.0113(16) 


48 


1 






0.2029(24) 


0.0391(21) 


0.0084(18) 


64 


1 






0.1946(31) 


0.0354(19) 


0.0072(19) 


96 


1 






0.1748(31) 


0.0256(26) 


0.0032(18) 


16 


2 


0.08486(17) 


0.11735(30) 


0.2187(13) 


0.0531(8) 


0.0140(7) 


24 


2 


0.08255(20) 


0.11505(32) 


0.2000(16) 


0.0436(10) 


0.0101(7) 


32 


2 


0.08146(20) 


0.11373(32) 


0.1893(13) 


0.0377(10) 


0.0088(8) 


48 


2 


0.08109(23) 


0.11337(37) 


0.1712(14) 


0.0306(10) 


0.0071(9) 


64 


2 


0.08051(23) 


0.11234(35) 


0.1629(18) 


0.0266(8) 


0.0026(9) 


96 


2 


0.08087(26) 


0.11337(44) 


0.1470(20) 


0.0207(11) 


0.0038(10) 


16 


4 


0.17102(20) 


0.22671(30) 


0.2088(6) 


0.0624(5) 


0.0240(4) 


24 


4 


0.16283(18) 


0.21812(28) 


0.1883(6) 


0.0479(5) 


0.0144(3) 


32 


4 


0.15965(18) 


0.21477(26) 


0.1755(7) 


0.0414(5) 


0.0107(4) 


48 


4 


0.15721(16) 


0.21206(23) 


0.1602(7) 


0.0338(5) 


0.0076(5) 


64 


4 


0.15602(20) 


0.21065(29) 


0.1489(8) 


0.0283(6) 


0.0071(5) 


96 


4 


0.15531(20) 


0.21012(30) 


0.1377(9) 


0.0234(6) 


0.0050(6) 



Table 7: Finite lattice renormalization group flow of the Aij for the dual of 
the XY model at fi = 1.1197 
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L 


/ 




A 2 ,i 


A 3J 


A 4J 


A 5 ,i 


12 


1 






0.2170(31) 


0.0474(22) 


0.0125(17) 


16 


1 






0.2061(23) 


0.0403(27) 


0.0068(15) 


24 


1 






0.1944(45) 


0.0407(32) 


0.0057(25) 


32 


1 






0.1729(41) 


0.0312(29) 


0.0058(27) 


12 


2 


0.08611(30) 


0.11901(45) 


0.1850(18) 


0.0409(12) 


0.0090(10) 


16 


2 


0.08408(27) 


0.11711(40) 


0.1725(16) 


0.0314(11) 


0.0062(12) 


24 


2 


0.08185(38) 


0.11418(58) 


0.1603(27) 


0.0275(14) 


0.0064(11) 


32 


2 


0.08140(34) 


0.11431(59) 


0.1458(27) 


0.0246(15) 


0.0051(13) 


12 


4 


0.17782(32) 


0.23433(45) 


0.1875(8) 


0.0588(6) 


0.0327(6) 


16 


4 


0.16839(24) 


0.22388(35) 


0.1663(8) 


0.0413(6) 


0.0143(6) 


24 


4 


0.16100(26) 


0.21601(40) 


0.1488(10) 


0.0305(8) 


0.0068(6) 


32 


4 


0.15800(25) 


0.21312(37) 


0.1380(10) 


0.0254(7) 


0.0050(5) 



Table 8: Finite lattice renormalization group flow of the Aij for the DG 
model at K DG = 0.6645 



L 


/ 


Ai.i 


A 2 ,i 


A 3 ,i 


A 4 ,i 


A 5 ,i 


32 


1 






0.2800(31) 


0.0766(27) 


0.0232(19) 


64 


1 






0.2404(35) 


0.0525(27) 


0.0132(25) 


128 


1 






0.2103(41) 


0.0425(28) 


0.0084(24) 


256 


1 






0.1866(45) 


0.0312(28) 


0.0028(24) 


32 


2 


0.08291(30) 


0.11513(52) 


0.2368(18) 


0.0609(12) 


0.0191(11) 


64 


2 


0.08202(33) 


0.11504(46) 


0.2023(25) 


0.0424(13) 


0.0088(12) 


128 


2 


0.08083(34) 


0.11337(52) 


0.1756(25) 


0.0334(14) 


0.0067(11) 


256 


2 


0.08069(33) 


0.11313(57) 


0.1553(26) 


0.0251(15) 


0.0035(13) 


32 


4 


0.16406(29) 


0.21936(45) 


0.2230(9) 


0.0656(7) 


0.0222(7) 


64 


4 


0.15914(27) 


0.21457(37) 


0.1891(9) 


0.0471(7) 


0.0136(7) 


128 


4 


0.15651(28) 


0.21137(43) 


0.1635(13) 


0.0320(7) 


0.0072(6) 


256 


4 


0.15513(28) 


0.21010(42) 


0.1412(11) 


0.0255(9) 


0.0057(7) 



Table 9: Finite lattice renormalization group flow of the Aij for the ASOS 
model at K ASOS = 0.8061 
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BCSOS 


model 




B 




^2,2 


B 




D 2 ,4 


6 


0.08313(15) 


0.11352(25) 


3 


0.17163(15) 


0.22080(22) 


8 


0.08244(21) 


0.11389(35) 


4 


0.16395(19) 


0.21711(29) 


12 


0.08175(22) 


0.11347(36) 


6 


0.15955(19) 


0.21389(29) 


16 


0.08106(23) 


0.11307(39) 


8 


0.15770(21) 


0.21210(31) 


24 


0.08078(22) 


0.11304(37) 


12 


0.15658(18) 


0.21132(28) 


32 


0.08077(17) 


0.11317(29) 


16 


0.15567(14) 


0.21054(21) 


48 


0.08067(20) 


0.11316(33) 


24 


0.15487(16) 


0.20959(25) 


64 


0.08033(21) 


0.11263(35) 


32 


0.15451(17) 


0.20908(27) 






Dual of XY mc 


Ddel 




B 




^2,2 


B 




D 2 ,4 


8 


0.08164(16) 


0.11379(29) 


4 


0.15854(19) 


0.21358(28) 


12 


0.08112(20) 


0.11347(32) 


6 


0.15725(17) 


0.21228(27) 


16 


0.08066(20) 


0.11285(32) 


8 


0.15651(18) 


0.21149(26) 


24 


0.08073(23) 


0.11298(37) 


12 


0.15581(16) 


0.21060(23) 


32 


0.08031(23) 


0.11212(35) 


16 


0.15523(20) 


0.20983(29) 


48 


0.08078(26) 


0.11327(44) 


24 


0.15496(20) 


0.20976(30) 


DG model 


B 




^2,2 


B 




D 2 ,4 


6 


0.08051(28) 


0.11275(43) 


3 


0.15632(28) 


0.21143(41) 


8 


0.08089(26) 


0.11356(39) 


4 


0.15610(22) 


0.21092(33) 


12 


0.08043(37) 


0.11262(57) 


6 


0.15548(25) 


0.21022(39) 


16 


0.08060(34) 


0.11342(59) 


8 


0.15489(25) 


0.20986(36) 






ASOS 


model 




B 




^2,2 


B 


D h4 


D 2 ,4 


16 


0.08209(30) 


0.11424(52) 


8 


0.16083(28) 


0.21601(44) 


32 


0.08182(33) 


0.11482(46) 


16 


0.15834(27) 


0.21374(37) 


64 


0.08078(34) 


0.11331(52) 


32 


0.15631(28) 


0.21116(43) 


128 


0.08068(33) 


0.11312(57) 


64 


0.15508(28) 


0.21005(42) 



Table 10: Results for the flow of the _D 8j ; at criticality 
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5.3 Determination of non-universal constants 

The matching also allows us to determine the non-universal constants appearing in the 
formulae for the divergence of observables near the roughening transition. Let us discuss 
this over the example of the correlation length £. Its critical behaviour is 



£ ~ A exp (Ck" 1/2 ) , k 



K-Kr, 
Kr 



(39) 



Let us consider matching on an RG trajectory in the phase with finite correlation length 

(smooth phase), close to the critical trajectory. Let us assume that the BCSOS block 

observables match the SOS block observables for sufficiently large B BGSOS ', B sos , with 
B sos = b sos B Bcsos_ Then 

^sos = b sos^Bcsos_ (40) 
Inserting eq. (39) in eq. (40) we get 

f os ~ b s m os A BCSOS exp [c BCSOS (k bcsos ) " 1/2 ) . (41) 

It is a general assumption of the renormalization group that couplings on the blocked 
system are smooth functions of the block size and the coupling on the fine lattice. We 
thus assume that b^ s is a smooth function of k sos , even at the roughening transition, 

b s m os (K) = b s m os (K R ) + O (k sos ) . (42) 
Furthermore, also K BCSOS j s a smooth function of k sos , 

K BCSOS = qK sos + O ((« 5 ° 5 ) 2 ) • (43) 
In the limit k sos — > we get 

^ A sos exp (C 505 ^ 505 )" 1 / 2 ) , (44) 

with 

j^SOS _ fySOS j^BCSOS 

c SOS = q -l/2 c BCSOS_ (45) 

With the results obtained above for b m and with A BCSOS = | we find 

A XY = 0.223(13), 
A DG = 0.078(5), 
a asos = 0.70(8). (46) 
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Figure 7: at criticality, plotted as a function of B 
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Figure 8: Ai 2 at criticality, plotted as a function of B 
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Figure 9: 1^2,2 a t criticality, plotted as a function of B 
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Figure 10: A 2) 2 at criticality, plotted as a function of B BCSOS = B sos /b 
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Figure 11: at criticality, plotted as a function of i? 
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Figure 12: Ai 4 at criticality, plotted as a function of B 
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Figure 13: a t criticality, plotted as a function of i? 
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Figure 14: A2 4 at criticality, plotted as a function of B 
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Figure 15: at criticality, plotted as a lunction ol B 



BCSOS _ nSOS II.SOS 



0.26 
0.24 
0.22 
0.20 h 
0.18 
0.16 h 
0.14 



— 1 1 1 1 - 



22 



_l 1 1 l_ 



m 



- 1 1 1 1 1 1 r - 

'bcsos.a32' ne— 1 

'xy.a32' h — 1 . 

'dg.a32' ra-n 

'asos.a32' hx— 1 



35 



10 



B BCSOS 



100 



Figure 16: A3 2 at criticality, plotted as a function of B 
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Figure 17: A^ : 4 at criticality, plotted as a lunction ol B 
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Figure 18: A41 at criticality, plotted as a function of B 
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Figure 19: A^2 a t criticality, plotted as a lunction ol B 
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Figure 20: A4 4 at criticality, plotted as a function of B 
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Figure 21: at criticality, plotted as a function of B 
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Figure 22: A5 2 a t criticality, plotted as a function of B 
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There remains to determine the constant q connecting k sos and k bcsos . For this purpose 
let us rewrite the matching condition for coupling constants in the neighbourhood of the 
critical point. In order to keep the formulae compact, we shall write an '5" for 'SOS' and 
a '_B' for 'BCSOS'. The matching condition reads 



A*, [B> , K> R + k')= A?, [B» , K% + k») . (47) 
A Taylor expansion around I(f> and K B , respectively, yields 

* M + * s {0i ( BS -^) + a S ( fiS - ^) S ("« ] 

= (same for BCSOS) . (48) 

These equations simplify since the observables match at criticality. As a consequence of 
the matching condition, 

A?,(fl s ,A-J) ^ Af,(fl*,A-| 

^(Ki) . 0. ,49) 

The second of these equations expresses the fact that we keep B s fixed when tuning the 
other parameters in order to fulfil the matching condition. We are left with 



w^v '^ 



/ \ 

• (50) 



dA ll ( nB „B\ . ( r>R „B\ 9B 



8KB \ BB ' K r) + Q^B \ BB ' K r) d^B^ 



Let us now again restrict our attention to i = 1,3. We can then solve eqs. (50) with 
respect to k s /k B : 



rl A B rl A B , , „ , , „ 

K _ dK B ' dB B dK B ' dB F 
k B dA? dA B 9A? dA B 



(51) 



_ I 1,1 _ 3,1 I 

dK s l dB B dK s l dB B 

The partial derivatives dA/dK can be determined with reweighting methods or with the 
help of the formula 

8 A 

= -(AH) + (A)(H) . (52) 

H is the Hamiltonian, and the expectation values are taken in the system with partition 
function Z = J2coni ex P( — & 'H). The partial derivatives dAfjdB B can be extracted from 
table 5. 

However, eq. (51) simplifies very much if the term dij in eq. (50) can be neglected. 
This is the case when 

fl uB 

:A'|)~0. (53) 



dB B 

~dK B ~ [ 



RG matching with the BCSOS model 



37 



Dual of XY model 


L 


i = 3, / = 1 


i = 1, / = 2 


i = 3, / = 2 


i = 1, / = 4 


i = 3, / = 4 


16 


0.413(8) 


0.430(11) 


0.406(8) 


0.424(10) 


0.400(7) 


24 


0.413(9) 


0.436(14) 


0.419(8) 


0.431(14) 


0.418(8) 


32 


0.427(12) 


0.461(19) 


0.429(10) 


0.437(19) 


0.421(9) 


48 


0.434(15) 


0.441(21) 


0.423(12) 


0.445(22) 


0.445(9) 


64 


0.449(17) 


0.420(28) 


0.445(13) 


0.447(29) 


0.434(11) 


96 


0.449(27) 


0.447(46) 


0.440(21) 


0.437(45) 


0.424(15) 


Discrete Gaussian model 


L 


i = 3, / = 1 


i = 1, / = 2 


i = 3, / = 2 


i = 1, / = 4 


i = 3, / = 4 


12 


0.381(12) 


0.394(17) 


0.375(10) 


0.394(18) 


0.369(9) 


16 


0.381(13) 


0.384(18) 


0.379(11) 


0.385(19) 


0.374(8) 


24 


0.388(17) 


0.384(30) 


0.381(13) 


0.388(29) 


0.374(11) 


32 


0.397(24) 


0.364(38) 


0.394(19) 


0.372(39) 


0.387(13) 


ASOS model 


L 


i = 3, / = 1 


i = 1, / = 2 


i = 3, / = 2 


i = 1, / = 4 


i = 3, / = 4 


32 


1.482(32) 


1.504(35) 


1.451(25) 


1.426(40) 


1.345(21) 


64 


1.385(42) 


1.474(46) 


1.345(40) 


1.457(44) 


1.447(42) 


128 


1.383(50) 


1.416(66) 


1.373(43) 


1.510(72) 


1.421(41) 


256 


1.909(106) 


1.625(157) 


1.712(75) 


1.608(151) 


1.586(57) 



Table 11: {k s /k B )^" e for the three SOS models 




Equivalently, the simplification relies on the assumption that the expansion of the matching 
B B around the roughening coupling K B is of second order in k B . Then we get 

If the approximation is a good one, the (k s fk B )^" e should be independent of Our 
results for these quantities are summarized in table 11. 

For all three models the numbers for (k s fk B )^" e for different are consistent 

within the error-bars. We arrive at the following estimates: 

k BCSOS /k XY = _ 43(1)j 
k BCSOS /k DG = _ 39(1)j 
k BCSOS /k ASOS = L46(6) _ (55) 
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Now we use that 

/ .,sos \ 1 / 2 



n K 

,SOS _ uSOS I is SOS 



r<sos _ ^ r ,Bcsos /r K \ 

° ~ 1 ...BCSOS ° 



and k = k bUb /Kf { Ub . We get 



C XY = 1.78(2), 
c dg = 2.44(3), 
C A505 = L14 ( 2 ). (57) 



5.4 Comparison with other numerical studies 

We compared our results with those obtained in other Monte Carlo studies. 

Janke and Nather [51] simulated the XY model with the Villain action in the vortex 
phase using Wolff's single cluster algorithm [60]. They measured correlation lengths up 
to £ = 140 on lattices up to L = 1200, at (3 ranging from (3 V = 0.590 up to (3 V = 0.675, 
where fi v = 0.5/ K DG . They fitted their results for £ to eq. (39). To check for systematical 
errors due to a too large distance to the critical point, they used two different definitions 
of k: 

kt = \T-T C \/T C , 

k p = \(3- (3 C \I(3 C . (58) 

The two definitions agree to the first-order Taylor expansion around the critical point. 
Hence both fits should give consistent results when the data included are obtained in a 
sufficiently small neighbourhood of the critical point. The comparison of the results of 
Janke and Nather with our results is given in table 12. The results of the two fits are not 
consistent within the error-bars. One is therefore led to the conclusion that the systematic 
error due to a too large distance of the simulation points from criticality is much larger 
than the quoted statistical errors. 

The authors give (3 C = 0.752(5) as an overall estimate for the critical coupling. Taking 
into account the systematic errors of the fits, the results of their simulation are well 
consistent with our results. 

In [30], the roughening coupling of the DG model was estimated from fits of the finite 
size behaviour of the surface thickness. The fits were done with a renormalization group 
improved formula. The best estimate was (3 C = 0.5/Kr = 0.755(3), which is nicely 
consistent with the estimate arrived at in the present paper, namely (3 C = 0.7524(7). 

In the case of the XY model with cosine action we can compare our results with a 
fit given in [75], which includes data of [50] and data of the authors, see table 13. As in 
the case of the DG model one can say that the results compare well with ours taking the 
systematic errors of the fits into account. We also include the results of the MCRG study 
[52] in this comparison. 



Solution of the KT equations 



0.06 
0.05 - 
0.04 - 
0.03 - 
0.02 - 
0.01 
0.00 



- 1 1 1 1 1 1 - 



m 



—i 1 



'bcsos.a54' ho — i 

'xy.a54' h — i 

'dg.a54' i-b— i 

'asos.a54' ^-h 



m 



_l I L_ 



10 

B BCSOS 



Figure 23: A5 4 at criticality, plotted as a function of B 
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B bUb lbl 



Definition of k 






C 




0.75106(36) 


0.1204(18) 


2.370(11) 




0.75814(40) 


0.0287(7) 


2.812(14) 


This work 


0.7524(7) 


0.078(5) 


2.44(3) 



Table 12: Comparison of our results for the DG model with those of ref. [51]. 
The relation between (3 V and K DG is (3 V = 0.5/ K DG 



Authors 




A 


C 


Gupta et al. 


1.1218 


0.2129 


1.7258 


Biferale 


1.112(2) 




1.74(20) 


This work 


1.1197(5) 


0.223(13) 


1.78(2) 



Table 13: Comparison of our results for the XY model with those of refs. [75] 
and [52] 
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Appendix 1: Solution of the KT equations 

In this appendix, we will solve the KT equations 

y =-xy, 

x =-y 2 . (59) 

It is easy to see that dy 2 jdt = dx 2 jdt. Therefore E = y 2 — x 2 is a constant of motion. We 
consider separately the three cases E = 0, E < 0, and E > 0. 

E = 0: From y 2 — x 2 = it follows that either x = y or x = —y. In the first case one 
is on the critical trajectory that runs into the fixed point at (x,y) = (0,0). The second 
case corresponds to the expanding manifold that leaves the fixed point. On the critical 
trajectory y obeys the differential equation y = —y 2 , which has the solution 

y(t) = XT7 ' (60) 

yo 

where yo = y(t = 0). On the expanding manifold y obeys y = y 2 , and the solution is 

f(*) = XT7- (61) 

yo 

Note that on this trajectory y diverges for finite /. 

E < 0: We look at the solutions in the x > region (the solutions for x < can be 
obtained analogously). Let a denote the point where the trajectory hits the s-axis at 
infinite /. Define x = a + / . We then have the system of differential equations 



dy 2 
dt 



-2(a+f)y 2 , 



f =-y 2 - (62) 



It follows that 



and therefore 



^ = 2(a + /), (63) 



y 2 = f 2 + 2a f . (64) 



The integration constant vanishes because y = for / = 0. We find the following diffe- 
rential equation for /: 

| = -/(/ + 2a). (65) 
This equation can be solved by separation of variables 

_J_V ,66) 
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By integration we find 

/W = 2a{(l + 7 ^)e 2 ^)-l}" 1 . (67) 
E > 0: Define e = \f~E~. From y 2 = e 2 + x 2 one obtains 

This differential equation can again be solved by separation of variables. The solution is 

x(t) = e tan ^arctan — — - — e(t — /q)^ • (69) 

Appendix 2: Gaussian model block correlation functions 

In this appendix, we shall derive an explicit formula for the block correlation functions of 
the massless Gaussian model defined by the Hamiltonian 

U{y) = \ (¥>, -Atp) = ^) 2 • (70) 

<%,y> 

The block lattice A' consists of / X / blocks x' . Each block x' is built of B X B fine grid 
points x, so that the fine lattice has extension L X L with L = IB. All lattices involved are 
provided with periodic boundary conditions. We define a block-block correlation function 

by 

G(x',y') = < {<t> x , - <t> y ,) 2 > , 

= B~ 2 J2^- (71) 

x£x' 

G has the following Fourier representation: 
G(x\y') = 2B-H- 2 Y J p- 2 IE E ex P(^-^))- E E ex P(^-y))j • (72) 

PT^O \x{zx ! y£x' x£x' y£y' / 

The periodic momentum squared is defined through 

ft 2 = 4 — 2 cospi — 2 cosp2 5 (73) 



where the pi denote components of the lattice momentum p. The pi run over the values 
L 



(0, 1, — 1) ^r. It is not difficult to derive that 



E E ex P (^( s ~ 2/)) = F M ^(^2) ex P (ipB(x' - y')) , (74) 

x£x' y£y' 
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with 



F( P i) 



B 2 if Pi = 0, 



(75) 




otherwise . 

\Pi) 



We thus arrive at 



G(0,x') = 2B~ 4 L 



2 Y J P~ 2 F(Pi)F(P2) (1 - cos(px'B)) . 



(76) 



The effective Laplacian defined in section 4 can be obtained as follows: 



A eff (0, x') = l~ 2 exp(iPx') A eff (P) 



(77) 



p 



Here the summation over all block lattice momenta Pi = (0, f, ...,/— and A e ff(P) is 

given by 
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